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ABSTRACT 

Vie present theoretical models for the Ionization of trace elements In a 
strong stellar wind by a compact binary X-ray source and for the resultlna 
orbital phase dependence of the emergent soft X-ray spectra and the profiles 
of ultraviolet resonance lines. Model results aqree qualitatively with the 
X-ray and ultraviolet spectra of the system 4U 0900-40/HD 77581 and exolaln 
the suppression of the absorption profiles of the SI IV A1394 and C IV U548 
lines when the X-ray source Is In front of the star. The model predicts that 
the absorption profiles of the N V *1239 and 0 VI *1032 lines will be enhanced 
rather than suppressed during this orbital phase. We predict phase-dependent 
linear polarization In the resonance line profiles. We discuss how future 
observations of these phase dependent effects In early-type binary X-ray 
systems may be used to Investigate the dynamics of stellar winds and their 
Interactions with the X-ray source. 

Subject Headings : X-rays: binaries - stars: winds - radiative transfer 


*0n leave from JILA and Department of Astrophysics, University of Colorado. 
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I. INTRODUCTION 

Early-type giant and supergiant stars have strong stellar Minds, with mass 
loss rates In the range M «• 10* 8 to * 10“® Mg yr'** and terminal velocities in 
the range « 500 to ■ 3000 km s' 1 (Conti and McCray 1980). The winds are 
probably driven by resonance scattering of the stellar radiation by Ions In 
the wind, but this Interaction Is not fully understood. In particular, the 
theoretical models do not accurately account for the Ionization of trace 
elements in the wind or for the details of the observed resonance line shapes 
and strengths, and they do not explain the observations of soft X-ray emission 
from such stars, which Is probably produced by shocks or instabilities In the 
wind (Lucy and White 1980). 

Several galactic 08 stars have compact X-ray source companions In close 
binary orbits. In such systems the winds can have sufficient column density 
to cause soft X-ray absorption variations with orbital phase (Buff and McCray 
1974). In turn, a luminous X-ray source can substantially modify the 
Ionization of trace elements In the wind, causing a phase dependence of the 
ultraviolet P- Cygnl type line profiles formed by these trace elements 
(Hatchett and McCray 1977). Observations have confirmed the existence of 
phase dependences in the X-ray and UV spectra from a few of the X-ray binaries 
with OB supergiant companions (e.g.., Dupree et al. 1978, 1980; Hammer schl ag- 
Hensberge et al. 1980; van der Klis et al. 1982; Kallman and White 1982; and 
references therein). These observations provide a probe of the density and 
Ionization structure In the winds of these systems. 

A compact binary X-ray source in a stellar wind Is not only a passive 
probe of the wind; it Is also an active probe of the physical mechanism for 
accelerating the wind. The X-ray emission can modify the dynamics of the wind 
directly through radiation pressure and Indirectly through changing the 
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Ionization of trace elements and, herce, the coupling of the stellar radiation 
to the wind (e.g., MacGregor and Vltello 1982; Friend and Castor 1982; 

Frans son and Fabian 1982). These effects are complex and poorly understood. 

In principle, we should be able to gain a satisfactory understanding of 
these systems by constructing theoretical models for the Ionization structure 
and dynamics of the winds and calculating the resulting orbital variations of 
their X-ray and UV spectra that may be compared to the observations In 
practice, this goal Is distant, because we do not even understand t!>: 
structure and dynamics of stellar winds of single OB stars, much less of 
binary X- ray sources, and the properties of the stellar winds are sensitive 
to many physical parameters that are not well ’..r.own. For single stars these 
Include the wind density (mass loss rate, velocity law) and ionization (which 
Is affected by the diffuse X-ray spectrum In the wind). For binary systems, 
additional uncertainties arise from Intrinsic variability of the X-ray source 
and from lack of knowledge of the soft X-ray spectrum (which strongly affects 
the Ionization of trace elements in the wind). 

In this paper we consider a more Idealized model. In which the dynamics of 
the stellar wind are not affected by the X-ray source, so that the wind's 
velocity and density distributions are assumed to be spherically symmetric and 
identical to what they would be In the absence of the X-ray source. However, 
ws do calculate quantitative models for the Ionization of the wind and for the 
variations of the X- ray and UV spectra resulting from the orbital motions of 
the X-ray source for the example of 4U 0900-40/H0 77581 (Vela X-l), the best 
observed such system. In Sec. II we describe a heuristic model for the 
structure and Ionization of a stellar wind around an isolated 0-star. In 5111 
we describe the Ionization structure of a stellar wind with an embedded X-ray 
source and the resulting orbital variations of the X-ray spectrum, and in §IV 
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we describe the variations In the UV resonance line profiles. In §V we 
compare these model results with present observations of 4U 0900-40/HD 77581, 
we discuss the limitations of the model, and we discuss what may be learned 
from further observations of this and other such systems. 


II. THE UNDISTURBED STELLAR WIND 

As shown by Olson and Castor (1981), the detailed calculation of the 
Ionization equilibrium of several Ions of several elements In the expandlnq 
winds of normal 0- and B-type stars Is quite complex. Besides direct 
photoionization from the stellar flux, diffuse radiation from the optically 
thick wind also contributes significantly to the photoionization rates. Trace 
amounts of Ions with large Ionization potentials, such as 0 VI , are produced 
by Auger Ionization from X-rays produced In the winds o* normal stars (nison 
1978, Casslnelll and Olson 1979). 

In the spirit of approximation adopted for this paper, we shall try to 
Infer the background stellar wind conditions from observed spectra rather than 
attempt to construct a detailed wind model, as was done by Olson and Castor 
(1981). The background Ion fractions and Ionization rates In the undisturbed 
stellar wind have been estimated In order to fit the theoretical profiles to 
the observed profiles according to the method suggested by Castor and Lamers 
(1979). The optical depth alor.* the radius, R, from the stellar center Is 
written as a function of the velocity: 

*■(«) * T(1 + y)( 1 - - »)\ m 

where w(R) * v(R)/Va, Is the radial wind velocity normalized t.o the terminal 
velocity v„, w Q » w(R*) Is the wind velocity at the stellar continuum 
photosphere, and T and r are adjustable parameters. We assume that the wind 
velocity Is spherically symmetric about the primary and Is given bv 


w(R) * 1 - 0.98R*/R, 


( 2 ) 


which corresponds to w 0 * 0.02 and p ■ 1 In the notation of Castor and 
lamer*. Eq. (2) gives good results for most 0-star P-Cyqnl line oroflles 
(cf. Germany et al 1981) . 

The profiles of SI IV and C IV In HO 77581 have been observed at several 
orbital phases (Dupree et al. 1980). Although the observed profiles do not 
agree In detail with toe theoretical profiles In the atlas of Castor and 
Lamers (1979), fits of toe profiles observed during X-ray eclipse (♦ « 0) 
provide an estimate of the distribution of scattering Ions In the wind of 
HO 77581 In the absence of X-ray Ionization. The parameters of these fits are 
given In Table 1. 

Unfortunately, observations of the N V and 0 VI resonance lines of 
HO 77581 are not available. Therefore, we assume that the undisturbed wind 
profiles are the s Me as those of toe star <c-0r1 (B0.5Ia), which Is very 
similar to HO 77581 but does not have a compact companion. Gathler, Lamers, 
and Snow (1981) have analyzed Copernicus spectra of <-0r1 and have fitted 
theoretical line profiles using equation (1). The parameters of their fits, 
which will be adopted as appropriate for HO 77581, are given In Table 1. 

In order to deduce the Ion fractions, the radial optical depth In eo. (1) 
must be equated to the Sobolev radial optical depth: 

t a (R) * STT - Va"i (3) 

where * 0 Is the rest wavelength of the transition, f A Is Its oscillator 
strength, and dv/dR Is the radial velocity gradient of the wind at the surface 
where a « * 0 (l-v/c). The density, nj, of Ions of element 1 In charge state 1, 
can be related to other quantities by using mass conservation in spherical 
symmetry: 
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(R) ■ n(R) * MA < g|/(4irR 2 uv(R)), (4) 

where n(R) Is the atomic density, Af Is the element abundance relative to 
hydrogen, Is the Ion fraction relative to the element, and u Is the mean 
molecular weight. Combining equations (1) through (4) will specify the 
Ionization fraction as a function of radius as determined from the fitted 
profiles. 

Because the contact X-ray source provides an additional source of 
photoionization. It Is convenient to determine an effective Ionization rate 
for tiie undisturbed wind. Balancing recombinations and Ionizations gives 

n 1 +l n e * "J *!• 

where n e and T e are the electron density and temperature, a (cm3 s-1) is the 
recombination rate coefficient, and C Is the effective Ionization rate. 

We assume that HO 77581 has a radius R* • 35R 0 , a wind temperature 
T e » 0.9 Teff * 23000 K, and a terminal wind velocity v. « 17nn km s-1. Since 
most of the carbon and silicon are In the forms C V and SI V, eauations fl-5) 
Imply that 

?(Si IV) » 0.073 mJ C(l-w) 2 x 4 w 2 dw/dx]’ 1 s" 1 (*) 

and C(C IV) - 0.12 M.gCU-w^xVdw/dxl- 1 s' 1 , (7) 

where M _g Is the mass loss rate In units of 10“ 6 M g yr' 1 and x ■ R/R* Is tHe 
distance from the primary In units of the radius of the stellar photosphere. 
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These "background" Ionization rates In the unperturbed stellar wind will he 
added to the rates produced by the compact X-ray source. 

We infer the Ionization fractions of N V and 0 VI from the line profiles 
of x-Orl, using equations (1-4), the stellar parameters and those In table 1. 
The results are: 

g(N V) » 0.011 fCg (1-0.98 /x) x' 1/2 (8) 


and 


g(0 VI) - 0.0041 (l-0.9fl/x)x‘ 1/2 . («) 

Auger Ionization produces all of the 0 VI seen In the winds of the 0- and R- 
type stars and the N V seen In early B-type stars (Casslnelll and filson 1979). 
Auger Ionization couples several Ionization stages, maklno It difficult to 
specify a meaningful effective Ionization rate. Therefore, since most of the 
nitrogen and oxygen are In the forms N IV and 0 IV, we take the results of 
equations (8) and (9) to represent the lower bounds of the abundances of v v 
and 0 VI that will occur In the wind. 
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III. IONIZATION AND X-RAY TRANSFER 

An accurate description of the Ionization and heat’-g produced by an X-ray 
source In a stellar wind Is complicated both by the physics of X-ray transfer 
and by the binary geometry. X rays have the property that they can completely 
Ionize all abundant trace elements, and that they can heat gas to temperatures 
greater than 10 7 K. The binary geometry requires a two- dimensional solution 
of the coupled transfer and Ionization- heating problem. In this Section we 
describe our technique for modeling the effects of a compact X-ray source 
embedded In a stellar wind and discuss some of the Important features of these 
models. 

An understanding of the Ionization structure of the X-ray Illuminated wind 
can be gained by considering a simplified model. If the effects of background 
Ionization due to the stellar radiation field and other sources are neglected 
and the wind Is assumed to be optically thin, then the Ionization state and 
temperature of the gas at any point depend only on the X-ray source spectrum 
and the Ionization parameter, ? * Lx/nr 2 , where L x Is the total X-ray 
luminosity, r Is the distance from the X-ray source, and n(r) Is the local gas 
density (Tarter, Tucker, and Sal peter, 1969). If the wind has a spherically 
symmetric density distribution, n(R) • R" 2 (steady flow at constant velocity), 
the surfaces of constant £, and hence of constant Ionization and temperature, 
are either spheres surrounding the X-ray source or open surfaces surrounding 
the primary (Hatchett and McCray, 197B). Generally, the geometry of these 
surfaces can be characterized by the dimensionless parameter 

q « en x D 2 /L x • R 2 v(R)/[r 2 v(0)], (10) 

where n x Is the wind density at the X-ray source and 0 Is the orbital 
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separation. The surface q » 1 determines the boundary between the open 
(q < 1) and closed (q > l) surfaces. Furthermore, Hatchett and McCray showed 
that these conclusion* are essentially unchanqed when X-ray absorption effects 
are considered. 

In the region of the wind where Ionization and heating are due both to 
X rays and to background Ionization by the stellar radiation field and other 
sources, an accurate description of the state of the gas requires the 
construction of more detailed models. In the absence of X rays the background 
ionization depends on R, the distance to the primary, and the dependence of 
the gas state on both X-rays and background Ionization breaks the simple i- 
scallng arguments described above. We have modeled the Ionization and 
temperature structure produced by a compact X-ray source In a stellar wind by 
solving the coupled Ionization equilibrium, thermal equilibrium, and radiative 
transfer equations along a series of rays from the X-ray source. The 
temperature. Ion abundances, and radiation field are computed as functions of 
distance along each ray uslnq the X-ray transfer code described by Kallman and 
McCray (1982). This code Includes the effects of a great number of atomic 
processes In a gas consisting of all Ion stages of the abundant elements H. 

He, C, N, 0, Ne, SI, S, and Fe. The radiative transfer equation Is solved In 
a simplified single-stream approximation. In which the diffuse emission at any 
point Is assumed to be produced by material with a spherically symmetric 
distribution centered on the X-ray source. This approximation will not lead 
to serious errors In these stellar wind models, since diffuse emission 1$ 
never Important In determining the Ionization state of the gas. We find that 
12 rays evenly distributed In polar angle with respect to the line of centers 
are sufficient to provide an accurate Interpolation of- the Ionization 
structure between the rays. 
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Input parameters for rite model Include the X-rav source luminosity and 
spectral shape. We use values appropriate to 4U 0900-40 (Becker et al . 

1978), consisting of a 20 keV bremsstrahlung spectrum with a total luminosity 
l x • io 36 erg s’ 1 . The 9 as density at any point In the wind is determined bv 
the total mass loss rate of M ■ 4 x 10“ fi vr"', based on examination of l<v 

line profiles during X-ray eclipse (Oupree et al. lQflO) and an assumed steady 
spherically symmetric flow with velocity law alven bv eouatlon (2) and 
v* • 1700 km s* 1 . The curves of constant Ionization parameter are shown In 
figure 1, where C ■ log^t. 

The abundances of Ions with threshold eneroles less than inn eV are 
affected by the copious UV flux from the prlmarv star. We model the ionization 
(hie to this consonant of the radiation by addlna to the flux at each point in 
the wind a blackbody field with temperature T* • 2fi0nn K, diluted accordlna to 
the distance from the primary. For SI IV and C IV the backoround Ionization 
rates (eqs. 6, 7) are added to the X-ray photolonlzatlon rates. For v v and 
0 VI we take the background Ionization fractions (eos. F, 9) as lower Sounds 
to the actual Ionization fractions when the X- rav flux Is weak. In order *n 
ensure that the Ionization structure >f other elements Far From the X-rav 
source approximates the true background Ionization In the wind we have imposed 
an artificial lower limit on the gas temperature of 90S of the hlackhodv 
temperature, l.e., T m f n * 23000 K (Klelfi and Castor 1978; fiathier et, al . 

1981). 

The Ionization structure of the stellar wind Is Illustrated in 'loures ?a-d, 
which show surfaces of constant Fractional Ion abundance, gj, For the Ions 
C IV, SI IV, N V, and 0 VI. Close to the X-rav source these surfaces closelv 
resemble spheres, as shown by Hatchett and McCray (197B), Tne X-ravs clearlv 
suppress the abundances of SI IV and C IV throuohout most of the hal*-sn*ce 


occupied by the X-ray source. Outside of this region other Ionization 
processes become Important, and on the far side of the primary the abundances 
are unaffected by the X-ray source. Since C IV and SI IV are less Ionized 
than the dominant Ions (C V and SI V), their abundances are suppressed near 
the X-ray source and Increase monotonlcally with distance from the X-ray 
source. Close to the primary star the high wind densities suppress the 
destruction of SI IV and C IV, causing their surfaces of constant abundance to 
partially wrap around the primary. 

In contrast, N V and 0 VI are more Ionized than the dominant Ions (N IV. 

0 IV, and 0 V), and their abundances (figures 2 c.d) are enhanced by more than 
two orders of magnitude throughout most of the half-space containing the X-ray 
source as a result of photoionization of the more abundant Ions. The 
enhancement of 0 VI Illustrates the penetrating power of the 1 keV X rays 
responsible for the production of this Ion via K-shell photoionization 
followed by Auger emission. Close to the X-ray source the abundances of N V 
and 0 VI are suppressed by Ionization to higher stages, and the surfaces of 
constant Ionization are nearly spherical. 

Our model calculations also yield the X-ray spectra emerging from the 
various rays from the source, displayed In figure 3. These spectra Illustrate 
the orbital variations that a distant observer will see as a result of the 
changing column density and Ionization state of the gas along the llne-of- 
slght to the source. Since the Inclination of the 4U 0900-40/HD 77581 binary 
system Is close to 90° (Conti 1978), we shall neglect the distinction between 
the angle 8 0 between the line of centers and the ray to the observer and the 
orbital phase angle and approximate 6 0 ■ *( 1 - 2 $), where orbital phase 4 ■ 0 
corresponds to the center of the X-ray eclipse. (The generalization to 
arbitrary orbital Inclination Is given by eg. A1 In Appendix A.) Hear X-ray 
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eclipse the line of sight to the X-ray source traverses a large column density 
of stellar Mind which Is not highly Ionized by X rays. The dominant opacity 
above 1 keV Is due to K-shell absorption by Ions of C, N, and 0, the cross- 
sections of which aro almost Independent of the number of L-$he11 electrons, 
and the absorption spectrum closely resembles that of neutral gas (e.g.. Brown 
and Gould 1970). For example, at orbital phase 4 » 0.2 the theoretical 
spectrum fits well with absorption by cold gas with column density 
H - 3.1 x 10 22 cm* 2 and agrees with the observed spectrum (Kali man and White 
1982 — but cf. SV). However, the X-ray opacity of the wind Is reduced 
significantly In the regions where C, N, and 0 are fully Ionized. This effect 
Is Illustrated by the dashed curves In figure 3, which show the emergent 
spectra that would result If the wind were neutral. For example, at t * 0*2 
the difference between the dashed curve and the solid curve corresponds to a 
reduction of column density of neutral gas by 40%. 

X rays can dramatically suppress the opacity at energies < 1 keV by 
photoionlzlng He II and all Ions of C, N, and 0 with L-shell electrons. As 
Implied by the results shown In figures 2 a- d, soft X rays can penetrate and 
ionize the lower gas densities at large distances from the primary, greatly 
reducing the soft X-ray opacity at the orbital phases farthest from eclipse. 
This effect Is apparent In figure 3 at $ * 0.5, where - 5% of the incident 
Intensity below 1 keV leaks out to large distances and the spectrum shows 
dramatic departures from the simple c d gas absorption (dashed curve). 
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IV. P-CYGNI LINE PROFILES 

In an Isolated 08 star, ultraviolet resonance lines are formed largely by 
resonant scattering of stellar continuum photons by Ions of trace elements In 
the stellar wind. A blueshlfted absorption component results from scattering 
of photons out of the beam from the stellar disk to the observer by Ions 
moving toward the observer. This absorption feature Is superposed on the 
stellar continuum plus an emission line resulting from scattering of Isotropic 
stellar continuum photons toward the observer by ions In the radially 
outflowing wind. Since the emission line Is formed by ions in the far side of 
the stellar wind as well as those In the near side, the emission component is 
roughly symmetrical about the resonance wavelength (except that some of the 
redshlfted emission is eclipsed by the primary). The superposition of 
absorption on the stellar continuum and emission results in a P-Cygni type 
profile, with a strong absorption feature on the blue side of the' resonance 
wavelength and a weaker emission feature on the red side. 

Asymmetries In the distribution of ions In the stellar wind resulting from 
the action of a binary X-ray source will manifest themselves as variations in 
the profiles of the P-Cygnl lines with orbital phase. One can see 
qualitatively (Hatchett and McCray 1977) that If an ion responsible for a 
given resonance line is removed by X-ray photoionization from a region of 
stellar wind on the near side of the star, the blueshlfted absorption will be 
suppressed. Conversely, If the abundance of an Ion is Increased on the near 
side, the absorption will be enhanced. 

Quantitatively, P-Cygnl line profiles can be calculated by using the 
Sobolev approximation (Sobolev 1957, 1960; Castor 1970; Lucy 1971; Castor and 
Laniers 1979), which takes advantage of the fact that the line shapes are 
formed by gradients In the velocity fields of these very supersonic winds 
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rather than by thermal Doppler motions. For a single resonance line, photons 
observed at a given wavelength are scattered only by gas within a narrow range 
of velocities, such that 

\ « x 0 a - v^c). (id 

where v, Is the component of wind velocity projected toward the observer. For 
monotonlcally accelerated winds the condition v, * constant defines a unique 
surface In the wind, called the resonant surface. Figure 4 shows these 
surfaces for a wind whose velocity profile Is given by eq. (2). 

We assume that the density and velocity structure of the wind are 
unaffected by the presence of the X-ray source, and that the Ionization of the 
wind can be calculated quasi statically, l.e., assuming balance between 
Ionization and recombination. As a result the wind ionization structure Is 
axl symmetric with respect to the line of centers of the binary system, and the 
line profile that will be seen by an external observer depends only on the 
angle, 8 0 , between the line of centers and the line of sight. 

The flux In a resonance line at a given wavelength Is calculated by 
integrating the specific Intensity emerging from such a surface over a family 
of rays parallel to the line of sight: 

L x * 4* / dn / m 0 pdp I x (p,n) (12) 

where p Is the Impact parameter of the ray with respect to the primary and n 
Is the azimuthal angle of the ray. The specific Intensity has the form 

I x (p,n) * exp(-T X )I c + [1 - exp(-t x )]S x , 


(13) 
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where the line optical depth, t x , and the source function, S x , are functions 
of the Ion density, the radial and transverse velocity gradients of the wind 
fit the resonant surface, and the angle, e, between the line of sight and the 
radial direction at the resonant point. Expressions for these functions and 
details of our method of evaluating eq. (12) are given In appendix A. Note 
that I x does not have azimuthal symmetry unless the observer Is viewing the 
binary system along Its line of centers, because the X-ray photoionization 
breaks the spherical symmetry of the Ion densities. 

I c Is the Intensity of the primary star's photosphere if the ray 
Intersects the primary; otherwise It Is zero. We have assumed that I c is a 
flat continuum. There is no obstacle to Including a photospherlc line profile 
or limb darkening (which may be significant — cf. Castor and Lamers 1979), 
but for simplicity we have not done so. We normalize our emergent line 
fluxes, l x , to the emergent continuum flux, L c * 4ir 2 R* 2 I c , to yield the line 
profile. 

Actually, the resonance lines considered are doublets, with the blue 
component having twice the oscillator strength of the red component. The 
doublet spaclngs are of the same order of magnitude as the Doppler shift 
correspcndl ng to typical stellar wind velocities. This circumstance 
substantially complicates the calculation of the line profiles, because it 
becomes possible for a photon to be scattered first In the blue component of 
the resonance doublet and again, at another location, by the red component. 

We have neglected this effect and have calculated the line profiles as if they 
were single lines with the parameters of the blue component. At the end of 
this Section we will discuss how we expect our results to be modified by the 
effects of the doublet separations. 

In order to understand the line profiles, it Is useful to refer to the 
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hatched curves In figures 2 a-d Indicating the surfaces where the radial 
optical depth, t x * 1 (eq. 3). Well Inside these surfaces the wind Is 
transparent to resonance scattering, negligible diffuse scattering of 
resonance lines Into the bean occurs, and few resonance line photons are 
scattered out of the bean. Since every point In the wind corresponds to a 
unique wavelength In the resonance line, one can Infer the wavelength range 
for which resonant scattering Is suppressed or enhanced by X-ray Ionization at 
phase + >0.5 by superposing figure 4 onto figures 2 a-d, and at other phases 
by rotating figure 4 through the corresponding angle and then superposing It 
on figures 2 a-d. 

For example, consider the phase dependence of the SI IV X1394 line for our 
model of the 4U 0900-40 /HO 77581 system. Illustrated In figure 5a. If there 
were no X-ray source In this system the theoretical SI IV line profile In this 
system would consist of the solid emission peak labeled <p « 0.5 on the red 
side of the resonance wavelength and the dashed absorption trough labeled 
t ■ 0, 0.25 on the blue side. Figure 2a shows that the zone of transparency to 
resonance scattering fills most of the half-space on the side of the X-ray 
source, and the superposition of figure 4 shows that this zone of transparency 
extends from v,/v. » -0.15 to V|/v„ » -1. Therefore, when the X-ray source Is 
on the near side of the star (<p >0.5), the absorption Is completely removed 
for sx/SXj < -0.15. The little peak at 5X/5Xy « 0.15 Is formed by forward 
scattering by SI IV Ions at large distances on the near side. When the X-ray 
source Is eclipsed ($ >0), SI IV Is again present In front of the star and 
the absorption trough returns, but the redshlfted wing of the emission Is 
suppressed slightly as a result of the photoionization of SI IV on the far 
side of the star. Note that tne absorption trough Is unaffected by the X-ray 
source even at quadrature (♦ « 0.25) because the transparent zone does not 


occult the star. 

The behavior of the C IV X1548 resonance line (figure 5b) Is similar, 
though less pronounced because the zone of transparency Is smaller. (C IV Is 
more durable than SI IV because of Its higher Ionization potential.) The 
surface for t x ■ 1 Is open, but the surface for t x * 1/2 (dashed curve) Is 
closed and roughly equal In size to the primary star. If there were no X-ray 
source the line would consist of the solid redshlfted emission profile and the 
dashed {$ » 0) absorption trough. When the X-ray source Is In front of the 
star (+ * 0.5). the zone of transparency occults the star and the absorption 
Is suppressed. This suppression Is greatest for 6X/5Xy ■ -0.53, corresponding 
to the Doppler shift where the zone of transparency has the greatest 
diameter. (Note that this blue-shift Is substantially greater than the value 
5X/5Xy » -0.33 correspondl ng to the wind velocity at the X-ray source.) 
However, when the orbital phase has changed by 30° ($ * 0.42) the zone of 
transparency occults only part of the primary, and by 60° ($ » 0.33) the 
absorption trough has returned completely. The zone of transparency Is not 
large enough to cause a noticeable suppression of the redshlfted emission at 
any phase. 

In contrast, the absorption profiles of the N V X1239 and the 0 VI X1032 
lines are enhanced at 4 « 0 r .5 as a result of the enhancement of the abundances 
of these ions by X-ray Ionization (figures 2c, d). In the absence of X rays (on 
the eclipsed side) the wind becomes optically thin (t x * 1) at R > 2.7 R* 

(v > 1100 km s“l) for N V X1239 (figure 2c) and at R > 9 R+ (v > 1500 km s“*) 
for 0 VI X1U32. For each case the resonance line profile in the absence of X- 
ray photoionization would consist of the solid U * 0.5) emission peak and the 
dashed ($ * 0) absorption trough. Note that this profile Is weaker for N V. 
owing to the lower abundance of nitrogen. Consequently, the effects of X-ray 


photoionization art wore pronounced for the N V line. The X-ravs 
substantially Increase the absorption at p • 0.5 and even at auadrature 
U » 0.25), while they enhance the emission at X-ray eclipse (♦ ■ 0). The 
suppression of the N V absorption due to the transparent zone near the X-ray 
source Is barely noticeable because this zone covers onlv a small fraction of 
the stellar disk. The effects of X-ray photoionization on the 0 VI line are 
similar, but less pronounced because the abundance of 0 VI without x-ray 
enhancement Is high enough to make the line profile almost saturated at all 
wavelengths. 

We now consider how these results should be modified to account for the 
doublet splitting of the resonance lines. The line sollttlnq and wavelengths 
of the lines are Indicated om figures 5a, d for each Ion, assumlna that 
v„ ■ 1700 km s" 1 . Crudely, one can guess that the actual line profiles will 
be some sort of superposition of two such profiles, the stronaer blue one 
shown centered at x^ (where * 0) and the weaker red one centered at x R ; 
but the superposition Is not linear. However, the profile of the blue 
component of the doublet will he unaffected by scatterlno In the red comnonent 
for x < x R - 6Xy. Thus, for example, In SI IV the predicted profile of the 
blue (1393.76 a) component will be unaffected for x < 1304.87 A, or 
6 a/6Xj < 0.14, However, for 0.14 < 6x/6Xy < 1 the red emission profile of the 
blue component will overlap with the blue absorption trouah of the red 
component, and we cannot accurately predict how the profile should behave In 
this range without a more detailed calculation. The situation is 
qualitatively different for the case of C IV, where the doublet soaclnq Is 
relatively small compared to dXy, and one can oredlct that the actual profile 
will behave roughly as a singlet, but smoothed out by ax « ?.50 A 
(Sa/oXj ■ 0.28). The worst case Is that of N V, where the doublet soaclna Is 
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comparable to the Doppler shift, (Xr - Xg)/5Xj * 0.57, and the only part of 
the profile which can be predicted with confidence Is the blue absorption wing 
for 6\/s\j < -0.43. The example of 0 VI Is similar to that of SI IV, In that 
the doublet spacing Is almost equal to the maximum Doppler shift, so that the 
predicted behavior of the absorption trough of the blue component should be 
unaffected by the line splitting. 


21 


V. DISCUSSION 

For binary X-ray systems with early-type companion stars we have predicted 
that both tha soft X-ray spectrum and the profiles of the ultraviolet 
resonance lines *111 vary with orbital phase as a result of photoahsorpt.lon of 
X-rays by Ions In the stellar winds of the companion stars. Some of these 
effects, such as the variable soft X-ray absorption and the suppression of the 
absorption troughs of the C IV U548 and SI IV A1394 resonance lines by the X- 
ray source, have been observed In a number of such systems, Includlna dij nonn-an, 
SMC X-l, LMC X-4, and Cyg X-l. Otter predictions, such as the enhancement of 
the absorption troughs of the N V A1239 and 0 VI A1032 resonance lines by the 
X-ray source, have not been observed yet. We no* consider what can be learned 
about these and otter such systems from comparison of the present theory with 
observations and from future Improvements In theory and observations. 

The theoretical models embody two main approximations. The first, that of 
neglecting the doublet separations In the calculations of the ultraviolet line 
profiles, Is not a serious problem In principle. It can be removed by a 
straightforward, although technically complicated, extension of the present 
calculations. The second approximation, that of neglectlno the effects of the 
X-ray source on the dynamics of the wind, presents a much more difficult and 
Important problem. As discussed In $ I, there are a number of wavs In which 
the X-ray source can affect the wind, and none are understood well. Since we 
are not confident of our ability to calculate these effects from first 
principles, we consider the present theoretical calculations to provide onlv a 
qualitative guide to the observations. However, we hope that the theory will 
provide some Insight Into what the observations tell us about the actual 
dynamics of the winds. 

We have modeled the best-observed such system, 4U noorMO/HD 775P1. we 
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have already mentioned In Sill that the theoretical soft X-ray absorption 
spectrum agrees well with the observed spectrum at orbital phase <i> » n.2, 
shortly after egress from eclipse. However, the observed soft X-rav spectrum 
does not va»y with orbital phase In the way predicted by the model. In 
particular, the soft X-ray spectrum observed before Ingress (♦ ■ 0.8) 
Indicates a much greater column density of absorbing matter than that at 
egress (Kallman and White 1982). This observation clearly suaoests some sort 
of density enhancement In the gas trailing the X-ray source. However, the 
structure of this density enhancement cannot be Inferred from the X-ray 
observations alone. It might be a small density enhancement or gas stream 
near the X-ray source or a much larger wake In the wind. More detailed X-ray 
spectroscopic observations, with better phase coverage, will a /low us to map 
the azimuthal distribution of column density In this structure, but thev will 
not tell us Its size. 

The ultraviolet resonance line profiles are formed by structures In the 
Ionization of the wind that are comparable In size to the companion star and 
the orbital separation. The absorption profiles of the SI IV and C IV are 
suppressed by ibout the amount predicted when 4U 0900-40 Is In front of 
HO 77581, Indicating Ionization zones In the wind of about the sl 2 e predicted 
by the theory. However, the profiles of the observed lines do not aoree In - 
detail with the theoretical ones. (This Is sometimes true even for the 
resonance line profiles formed In the winds of single stars.) Some of the 
discrepancy might be attributed to our neglect of the doublet effect In the 
theory. For example, In C IV a 1548 the observed absorption trough at <p * 0.5 
Is double, as predicted by theory, but the deeper absorption occurs at low 
velocity rather than at high velocity. This might be explained by the fact 
that both the 1548 A and the 1551 A components can contribute to the low 


velocity feature, while only the 1548 A component can contribute to the hlqh 
velocity feature. Some of the discrepancy nrfqht also arise from an unknown 
photospherlc contribution to the resonance line. 

Observations of the resonance absorption line profiles offer the 
opportunity to Infer the site of the density structure In the wind that Is 
responsible for the asymmetries In the orbital variations of the soft X-ray 
spectra. If this structure Is comparable In size to HO 77581, asymmetries 
should also show up In the orbital variations of the ultraviolet lines; but If 
It Is much smaller than HD 77581, It should not significantly affect the 
ultraviolet lines. It Is possible, but not at all certain, that the large- 
scale density structure of the stellar wind Is not affected much by the X-rav 
source. Abbott (1981) has suggested that the radiation force on the wind due 
to resonance scattering of stellar continuum photons should be fairly 
Independent of the Ionization state of the wind, because so many different 
Ions and resonance lines can contribute to this force. If so, the dynamics of 
the stellar wind would be modified by X-ray photoionization and qravlty only 
In a region very near the X-ray source, much smaller than the comoanlon star. 

We do not expect the red emission components of the P-Cygnl lines to 
provide very useful diagnostics of the X-ray Influence on the wind, for three 
reasons. These are: first, the emission component Is built of small 
contributions from a region much larger than the 0 star, and therefore the 
effect of changing the Ionization state In a moderate volume near the n star 
Is small. Our calculations Illustrate this. Second, we know from observations 
of single 0 stars that the emission components do not agree well with rhe 
profiles calculated by the Sobolev model In which a flat continuum Is assumed 
for the 0 star's photosphere. The observed emission strength Is often much 
less than that from the Sobolev calculation. This Is not well understood; 


various possible explanations have been discussed by Abbott, Bohlln and Savaoe 
(1982), of *h1ch the most likely Is that the photospherlc continuum Is In 
reality a broad absorption feature. Third, the emission component Is formed 
by scattering In both components of the doublet. 

Although observations have demonstrated the possibility of uslnn 
observations of contact X-ray binaries to probe stellar *lnds, the oresent 
data have serious limitations. In particular, the X-ray sources are kno*n to 
vary Intrinsically by large (> 2) factors In timescales less than the orbital 
period. The resonance line profiles *111 vary as a result of these Intrinsic 
X-ray luminosity variations. Therefore, simultaneous observations of both the 
X-ray spectra and the ultraviolet spectra as a function of orbital phase *111 
be required In order to distinguish the effects of Intrinsic X-ray variability 
from those of the changing orbital geometry. In addition, the present 
ultraviolet spectra, even for 4U 0900-40/HD 77581, the brightest such system, 
have Inadequate phase coverage and slgnal/nolse to permit a detailed analysis 
of the Ionization structure In the *1nd. 

Finally, *e comment on ho* the effects described here should deoend on the 
parameters of other binary X-ray systems *1th early-tyoe companion stars. The 
degree of Ionization of trace elements In the *1nd Is controlled larqelv by 

e 

the parameter S 0 » 4*uv(D)L x /M. Thus, systems *1th more luminous X-rav sources 
or *eaker *1nds than 4U 0900-40/HD 77581 *111 have relatively laroer 
Ionization zones and their ultraviolet resonance lines should vary more 
dramatically. For example, In such systems the Ionization zones of M V and 
0 VI might be large enough to cause their resonance lines to vary In the same 
*ay as the SI IV and C IV lines do. 

Similar orbital variations have been observed In the ultraviolet resonance 
lines from the early-type binary X-ray systems Cvo X-l (Treves etai. I°8n), 



SMC X-l ( Hammer schl ag-Hensberge et al. 1983), and LMC X-4 (van Her *Hs et al . 


1982). These system do seem to have more luminous X-rav sources and weaker 
winds than 4U 0900-40/HD 77581, and their resonance lines behave rouohiy as 
expected. Unfortunately, these system are so faint that the JUP spacecraft 
cannot obtain high resolution ultraviolet spectra with oood signal /noise and 
phase coverage. Detailed theoretical models do not seam warranted by the 
presently available observations. 

Good high dispersion ultraviolet spectra have been obtained with lijr for 
the early-type binary system 4U 1700-37/Hf) 153919 (Ouoree et ai. 1 07 ft 1 ; they 
show no discernible variations of the SI IV, C IV, and M V resonance lines. 
This result can be understood because the X-ray luminosity of this system is 
low and the stellar wind is strong, so that £ 0 Is small and the rones where 
the resonance line opacities are affected by X-ray photoionization are much 
smaller than the companion star. 

Another possible way to observe asymmetries In the distribution of Ions in 
stellar winds Is to measure linear polarization In the resonance line 
profiles. One can see qualitatively that the scattered radiation In the line 
should be linearly polarized with the electric vector parallel to the 
projected orbital axis when the X-ray source Is displaced to either side of 
the companion star. Since the production of polarized nhotons hv resonant 
scattering reaches a maximum when t "1 (Goldrelcb and Kvlafls 1992), the 
polarization should show up best at those wavelenoths In the line profile 
corresponding to projected wind velocities (fioure 4) in the regions where 
t *l (figures 2a-d). Thus, we predict that the maximum In the linear 
polarization will appear first In the redshlfted part of the line profile as 
(possibly before) the X-ray source emerges from eclipse and win move to 
shorter wavelengths, reaching Its maximum degree for «t 9 « 0 . 2 *, then 
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tht blueshlfted part of the profile. The polarization should vanish at * ■ 
0.5, and the sequence should repeat In reverse as the source moves throuqh $ - 
0.75. Although a detailed calculation Is beyond the scope of this paper, #e 
estimate that the linear polarization In the line profile may reach a fe« 
percent under optimal circumstances. 

It should be possible to observe these and other such systems In much 
better detail #1th future Instrumentation. The Space Telescope *in he able to 
obtain excellent high dispersion ultraviolet line profiles and measure the 
polarization #1thfn them, and the Advanced X- ray Astronomy Facility *111 be 
able to provide complementary X-ray spectroscopic observations. Such 
observations #111 be a powerful tool for exploring the dynamics of stellar 
rflnds and their Interactions #1th compact X-ray sources. 

This #o He *as supported In part by grants from NASA (NSG-712*) and NSF 
{AST 80-19874) to the University of Colorado and by computlno facilities 
provided by the National Center for Atmospheric Research, #hlch Is suoported 
by the NSF. 
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TABLE 1 

PROFILE PARAMETERS FOR UV RESONANCE LINES 


Line 

*(*) 

_T 

2 

it 

C IV 

1548.19 

20 

2 

3.3 (-A) 

N V 

1238.81 

1 

0.5 

9.1 (-5) 

0 VI 

1031.93 

2 

0.5 

6.6 (-4) 

Si IV 

1393.76 * 

6 

2 

3.3 f-q 



APPENOIX A: CALCULATION OF LINE PROFILES 


The geometry used for our calculation Is depicted In Flqure 6, where: R Is 
the radial distance from the center of the primary (0) to the point P on the 
resonant surface; r Is the radial distance from the secondary (X) to P; p Is 
the Impact parameter of the ray from P to Earth with respect to the ray From 0 
to Earth; n Is the azimuthal angle of the ray, is the anole between the 
line of centers 0X_and the line of sight from 0 to Earth; a Is the the anole 
between the radius vector to the resonant point OP and the line of siaht from 
P to Earth; Is the angle between 0P_ and OX; and b x Is the anole between nx_ 
and PX . We define the corresponding direction cosines: u * cosQ; u 0 * cnsQ 0 ; 
Mp * cos Op ; and u x ■ coso x . Fon arbitrary orbital inclination anole, 1, we 
have the relation 


cos o * - sin 1 cos 2*<p. 
o 

The line source function for the conservative scatterlna process that we 
assume Is determined by the rate of scatterlna of continuum Dhotons at that 
point, i.e.. 


s - i c e c /e, 


(A?) 


in which the escape probability functions S and t$ c (Castor 1°7D1 are oiven h v 
angle Integrals Involving the Sobolev optical depth for each direction at the 
resonant point: 


1 

W 


1-e 


-X 


du, and * 


1 r 1-e' 


du 


( A3 1 


T 
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The Integral for » c Is extended only over those directions that intersect the 
primary's photosphere. 

The Sobolev optical depth, t(u) is zero If the primary lies between the 
observer and the resonant point alonq the ray. Otherwise It depends on the 
llne-of-slght velocity gradient, which In turn Is a function of the local 
radial and transverse velocity gradients, u, and the Ion density at the 
resonant point. We have the relation 

t * t 0 /{1+om2) $ ( A41 


In which x Q is the transverse optical depth and o Is the velocity-gradient 
anisotropy factor defined by 


o 


R dv . 
v W " 


(AS) 


The transverse optical depth t q can be expressed In terms of the radial 
optical depth t rad by 


T o * T rad (1+a) - 

For the P-Cygnl lines that we are conslderlna the ooacltv is verv nearly 
proportional to the total number of Ions In question, since almost all such 
Ions are In the ground state. If we deslanate by a* the relevant Ionization 
fraction and If, furthermore, q^° Is the corresponding Ionization fraction In 
the absence of the X-ray source, then we can write 


(A7) 


Vad * <9i/9i 0 ' Vad°* 


rfhere * ra d°* the radial optical depth In the absence of the X-rav source, is 
given by eq. (1) and the parameters listed In Table 1. 

Our computational technique Is the following: equations (12) and (13) are 
put In the form 

Ly/Lg * 1 + y j (l-y 2 r 1/2 dy (j* (S-J.)(l-e -t )2pdp * /*2p 3 S(l-e" T )du}, 

' l (A.) 

In rfhlch *e have defined y ■ cos n and u » 1/n. Note that o 3 S (l-e' T ) tends 

to a finite limit as p becomes large and the freauencv Is held fixed. The 

Integrals In eq. (A8) are done using appropriate variations of Gaussian 
quadrature. There Is reflection symmetry *1th respect to the niane of the 
line of sight and the line of centers, so the Integration over n need onlv 
cover the range 0 to » (y ■ -1 to 1). For this Integral use Gauss- 
Chebyshev quadrature *1th 12 abscissae y^. (This Is the same as a trapezoidal 
rule versus n In rfhlch use Is made of the periodicity.) the Integral over p 

from 0 to 1 Is done *1th a 12 point Gauss-Jacobl quadrature, that Is, the 

factor p In the Integrand Is Included as a welaht function. The u inteorai is 
done rflth standard 12 point Gaussian quadrature, transformed from the Integral 
(-1,1) to the Interval (0,1). The points and weights for these schemes are 
given by Abramowltz and Stegun (1964) In § 26.4. 

Since the Ionization has axial symmetry about the line of centers, 
depends on R and n p . The Interpolation of q^/o^° at the resonant noint 
provides * ra( j through eq. (A7), and therefore also t q throuah en. 1AF). The 
functions e and defined by eq. (A3) depend on t q , o, and u c a (1-1/R^)^*, 

the cosine of the cone angle subtended by the primary. We calculate these 
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using an approximation suggested by Ryblckl (private communication). The 
source function S follows from eq. (A2), the optical depth t alone the line of 
sight Is given by eq. (A4), and these are all the pieces needed to evaluate 
the Integrals In eq. (A8). 

The calculation proceeds as follows. Me define the dimensionless 
frequency variable. 


so that o « +1 is the violet edge of the absorption trouoh in the *-Cvon1 
profile and a ■ -1 Is the red edge of the emission component. Me select 
values for * and (we choose multiples of 0.05 and 30°, respectively). For 
given ° and p we find the radius of the resonant point from 


« « w(R)[l-p 2 /R 2 )] 1/2 (A10) 

(radial dimensions are normalized to R* and the terminal velocity has been set 
to unity). With the chosen velocity law (eo. 2), a non-linear eouat.lon for 9 
results, which Is solved uslnq the Mewton-Raphson method. Havlno solved eo. 
(A10) for the radius, the direction cosine u follows from 


w * a/w(R) 


(All) 


The direction cosine u p then follows from the spherical law of cosines: 

M p * V + y. 


(A12) 
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The coordinates R and Mp locate the resonant point p with respect to t,h<? 
centers of the primary and secondary stars. Specifically, the polar 
coordinates of P with respect to the secondary (X) are 

r ■ (R 2 ♦ D 2 - 2RDu ) 1/2 U13) 

and 

m x « (D - MpR)/r. ( A14) 
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FIGURE CAPTIONS 


Fig. 1 - 
Fig. 2a — 


Fig. 2b ~ 

Fig. 2c - 
Fig. 2d — 
Fig. 3 — 

Fig. 4 ~ 
Fig. 5a — 


Fig. 5b — 
Fig. Sc - 
Fig. 5d — 
Fig. 6 - 


Contours of constant Ionization parameter, 5 » L -/nr 2 *or the 
4U 0900-40/HD 77581 system. The curves are labeled by 
<. * logics* 

Contours of constant fractional abundance, g^, *or SI IV. The 
curves are labeled by (-logiQqj*). curve with hatch marks 
Indicates the boundary where the resonance line opacity t x « i, 
separating the optically thick (hatched side) *rom the oot.lcaiiv 
thin (solid side) zones. 

Same as Fig. 2a but for C IV. The dashed curve 
represents 1/2 . 

Same as Fig. 2a but for N V. 

Same as Fig. 2a but for 0 VI. 

Emergent X-ray continuum spectra as ^unctions of orbital phase, 
i*. Dashed curves: spectra that would result If the wind were 

neutral . 

Contours of constant projected stellar wind velocity, v #> f 0r 
the 4U 0900-40/ HD 77581 system. The curves are labeled bv the 
ratio v |/v« « m/^ t . 

Predicted profile of the SI IV A1394 line as a function of 
orbital phase, », for the 4U (190O-40/H0 77581 svstem. Abscissa 
Is line shift normalized to the maximum oosslble Poooler s h 1*t 
■ v«/c. Wavelengths and Ooopler shifts of resonance doublets 
are Indicated, assumlna that the resonance wavelenath js that of 
the stronger (blue) component and that v„ » 17no km s* 1 . 

Same as Fig. 5a but for C IV A1548. 

Same as Flo. 5a but for N V A1239. 

Same as Fig. 5a but for 0 VI A1032. 

Binary system geometry. 0 Is the location of the center ot the 
primary, X 1$ the center of the secondary, and p Is the location 
of a parcel of stellar wind. The line of sloht to n is not 
necessarily In the y,z plane. 
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